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\£) ' We compute the quantum Langevin equation (or quantum stochastic differential equation) repre- 

senting the action of a quantum heat bath at thermal equilibrium on a simple quantum system. These 
equations are obtained by taking the continuous limit of the Hamiltonian description for repeated 
quantum interactions with a sequence of photons at a given density matrix state. In particular we spe- 
cialise these equations to the case of thermal equilibrium states. In the process, new quantum noises 
are appearing: thermal quantum noises. We discuss the mathematical properties of these thermal 
quantum noises. We compute the Lindblad generator associated with the action of the heat bath on 
the small system. We exhibit the typical Lindblad generator that provides thermalization of a given 
quantum system. 



I. Introduction 

The aim of Quantum Open System theory (in mathematics as well as in 
physics) is to study the interaction of simple quantum systems interacting with 
very large ones (with infinite degrees of freedom). In general the properties that 
one is seeking are to exhibit the dissipation of the small system in favor of the 
large one, to identify when this interaction gives rise to a return to equilibrium or 
a thermalization of the small system. 

There are in general two ways of studying those system, which usually repre- 
sent distinct groups of researchers (in mathematics as well as in physics) . 

The first approach is Hamiltonian. The complete quantum system formed by 
the small system and the reservoir is studied through a Hamiltonian describing 
the free evolution of each component and the interaction part. The associated 
unitary group gives rise to a group of *-endomorphisms of a certain von Neumann 
algebra of observables. Together with a state for the whole system, this constitutes 
a quantum dynamical system. The aim is then to study the ergodic properties of 
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that quantum dynamical system. This can be performed via the spectral study of 
a particular generator of the dynamical system: the standard Liouvillian. This is 
the only generator of the quantum dynamical system which stabilizes the self-dual 
cone of the associated Tomita-Takesaki modular theory. It has the property to 
encode in its spectrum the ergodic behavior of the quantum dynamical system. 
Very satisfactory recent results in that direction were obtained by Jaksic and Pillet 
([JP1], [JP2] and [JP3]) who rigorously proved the return to equilibrium for Pauli- 
Fierz systems, using these techniques. 

The second approach is Markovian. In this approach one gives up the idea of 
modelizing the reservoir and concentrates on the effective dynamics of the small 
system. This evolution is supposed to be described by a semigroup of completely 
positive maps. These semigroups are well-known and, under some conditions, 
admit a generator which is of Lindblad form: 

C(X) = i[H, X] + ^ J2( 2L i XL * - L *i L - X ~ XL*Li). 

i 

The first order part of C represents the usual quantum dynamic part, while the 
second order part of C carries the dissipation. This form has to be compared with 
the general form, in classical Markov process theory, of a Feller diffusion generator: 
a first order differential part which carries the classical dynamics and a second 
order differential part which represents the diffusion. For classical diffusion, such 
a semigroup can be realized as resulting of a stochastic differential equation. That 
is, a perturbation of an ordinary differential equation by classical noise terms such 
as a Brownian motion usually. In our quantum context, one can add to the small 
system an adequate Fock space which carries quantum noises and show that the 
effective dynamics we have started with is resulting of a unitary evolution on the 
coupled system, driven by a quantum Langevin equation. That is, a perturbation 
of a Schrodinger-type equation by quantum noise terms. 

Whatever the approach is, the study of the action of quantum thermal baths is 
of major importance and has many applications. In the Hamiltonian approach, the 
model for such a bath is very well-known since Araki- Woods' work ([A-W]). But 
in the Markovian context, it was not so clear what the correct quantum Langevin 
equation should be to account for the action of a thermal bath. Some equations 
have been proposed, in particular by Lindsay and Maassen ([L-M]). But no true 
physical justification of them has ever been given. Besides, it is not so clear what 
a "correct" equation should mean? 

A recent work of Attal and Pautrat ([API]) is a good candidate to answer 
that problem. Indeed, consider the setup of a quantum system (such as an atom) 
having repeated interactions, for a short duration r, with elements of a sequence 
of identical quantum systems (such as a sequence of photons). The Hamiltonian 
evolution of such a dynamics can be easily described. It is shown in [API] that in 
the continuous limit (r — > 0), this Hamiltonian evolution spontaneously converges 
to a quantum Langevin equation. The coefficient of the equation being easily 
computable in terms of the original Hamiltonian. This work has two interesting 
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consequences: 

- It justifies the Langevin-type equations for they are obtained without any 
probabilistic assumption, directly from a Hamiltonian evolution; 

- It is an effective theorem in the sense that, starting with a naive model for 
a quantum field (a sequence of photons interacting one after the other with the 
small system), one obtains explicit quantum Langevin equations which meet all 
the usual models of the litterature. 

It seems thus natural to apply this approach in order to derive the correct 
quantum Langevin equations for a quantum heat bath. This is the aim of this 
article. 

We consider a simple quantum system in interaction with a toy model for 
a heat bath. The toy model consists in a chain of independent photons, each of 
which in the thermal Gibbs state at inverse temperature /?, which are interacting 
one after the other with the small system. Passing to the continuous interaction 
limit, one should obtain the correct Langevin equation. 

One difficulty here is that in [API], the state of each photon needed to be 
a pure state (this choice is crucial in their construction). This is clearly not the 
case for a Gibbs state. We solve this problem by taking the G.N.S. (or cyclic) 
representation associated to that state. If the state space of one (simplified) photon 
was taken to be n-dimensional, then taking the G.N.S. representation brings us 
into a n 2 -dimensional space. This may seem far too big and give the impression 
we will need too many quantum noises in our model. But we show that, in all 
cases, only 2n chanels of noise resist to the passage to the limit and that they can 
be naturally coupled two by two to give rise to n "thermal quantum noises" . The 
Langevin equation then remains driven by n noises (which was to be expected!) 
and the noises are shown to be actually Araki- Woods representations of the usual 
quantum noises. Furthermore, the Langevin equation we obtain is very similar to 
the model given in [L-M]. 

Altogether this confirms we have identified the correct Langevin equation 
modelizing the action of a quantum heat bath. 

An important point to notice is that our construction does not actually use 
the fact that the state is a Gibbs-like state, it is valid for any density matrix. 

This article is organized as follows. In section II we present the toy model for 
the bath and the Hamiltonian description of the repeated interaction procedure. 
In section III we present the Fock space, its quantum noises, its approximation by 
the toy model and the main result of [API]. In section IV we detail the G.N.S. 
representation of the bath and compute the unitary operator, associated with the 
total Hamiltonian, in that representation. In section V, applying the continuous 
limit procedure we derive the limit quantum langevin equation. In the process, we 
identify particular quantum noises that are naturally appearing and baptize them 
"thermal quantum noises", in the case of a heat bath. The properties of those 
thermal quantum noises are studied in section VI; in particular we justify their 
name. In section VII, tracing out the noise, we compute the Lindblad generator of 
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the induced semigroup on the small system. In section VII, being given any finite 
dimensional quantum system with its Hamiltonian, we show how to construct a 
Lindblad generator, representing some interaction with a heat bath, such that the 
quantum system thermalizes. 

II. The toy model 

We describe here the physical model of repeated interactions with the bath 
toy model. 

The quantum system (we shall often call "small system" ) to be put in contact 
with the bath is represented by a separable Hilbert space Hs , as state space, and 
a self-adjoint operator H$, as Hamiltonian. 

The toy model for the heat bath is the chain 

(g) C N+1 

of copies of € N+1 , where N > 1 is a fixed integer. Each copy of C N+1 represents 
the (simplified) state space of a photon. By this countable tensor product we 
mean the following. We consider a fixed orthonormal basis {eo, ei, . . . , ejv} of 
€ N+1 , corresponding to the eigenstates of the photon (e being the ground state); 
the countable tensor product is taken with respect to the ground state eo- Together 
with this structure we consider the associated basic matrices a*-, i,j = 1, . . . , N, 
acting on € N+1 by 

a) e k = S ik ej 

and their natural ampliations to ®keiN* < C- N+1 given by 

a if]A — / a% j on t ne k-th copy of € N+1 
3 \ I on the other copies. 

The Hamiltonian of one photon is the operator 

N 
i=0 

where the 7^'s are real numbers. Here notice two points. 

We have assumed the Hamiltonian Hr to be diagonal in the chosen basis. 
This is of course not actually a true restriction, for one can always choose such 
a basis. Note that H R describe the total energy of a single photon, not the whole 
field of photon. For this we differ from the model studied in [AJ1]. 

The second point is that 70 is the ground state eigenvalue, it should then be 
smaller than the other 7$. One usually assumes that it is equal to 0, but this is not 
actually necessary in our case, we thus do not specify its value. The only hypothesis 
we shall make here is that 70 < 7$, for alii = 1, . . . N. This hypothesis means 
that the ground eigenspace is simple, it is not actually a necessary assumption, it 
only simplifies our discussion. At the end of section V we discuss what changes if 
we leave out this hypothesis. 
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Finally, notice that the other eigenvalues 7^ need not be simple in our discus- 
sion. 



When the system and a photon are interacting, we consider the state space 
Hs <8> & N+1 together with the interaction hamiltonian 

N 

H I = J2{Vi®a° i +V*®ai), 
i=i 

where the V^s are bounded operators on Hs- This is a usual dipole-type interaction 
Hamiltonian. The total Hamiltonian for the small system and one photon is thus 

JV 

H = H s I + I <g> H R + (Vi ® a° + Vi <g> 4) . 

i=i 

Finally, the state of each photon is fixed to be given by a density matrix p which 
is a function of Hr. We have in mind the usual thermal Gibbs state at inverse 
temperature (3 : 

where Z = tr(e~ ,3HR ), but our construction applies to more general states p. 

Note that pp is also diagonal in our orthonormal basis. Its diagonal elements 
are denoted by {/3q, Pi, ... , P n }- 

We shall now describe the repeated interactions of the system Tig with the 
chain of photons. The system Hs is first in contact with the first photon only 
and they interact together according to the above Hamiltonian H. This lasts for 
a time length r. The system Hs then stops interacting with the first photon and 
starts interacting with the second photon only. This second interaction is directed 
by the same Hamiltonian H on the corresponding spaces and it lasts for the same 
duration r, and so on... This is mathematically described as follows. 

On the space Hs <8> € N+1 , consider the unitary operator representing the 
coupled evolution during the time interval [0, r]: 

U = e~ lTH . 

This single interaction is therefore described in the Schrodinger picture by 

p ^UpU* 

and in the Heisenberg picture by 

X ^ u*xu. 

After this first interaction, we repeat it but coupling the same Hs with a new copy 
of (D N+1 . This means that this new copy was kept isolated until then; similarly 
the previously considered copy of (D N+1 will remain isolated for the rest of the 
experience. 

The sequence of interactions can be described in the following way: the state 
space for the whole system is 

]N* 
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Consider the unitary operator U k which acts as U on the tensor product of Hs 
and the k-th copy of <D N+1 , and which acts as the identity on all the other copies 
ofC^ 1 . 

The effect of the k-th interaction in the Schrodinger picture is 

p i-> U k pUl, 

for every density matrix p on Hs ®%j € n+1 . In particular the effect of the k first 
interactions is 

where V k = U k U k -i . . . U±. 

Such a Hamiltonian description of the repeated interaction procedure has no 
chance to give any non-trivial limit in the continuous limit (r — > 0) without asking 
a certain renormalization of the interaction. This renormalization can be thought 
of as making the Hamiltonian depend on r, or can be also seen as renormalizing 
the field operators a°, a l of the photons. As is shown is [API] (see the detailed 
discussion in section III), for our repeated interaction model to give rise to a 
Langevin equation in the limit, we need the interaction part of the Hamiltonian 
to be affected by a weight 1 / \fr. Hence, from now on, the total Hamiltonians we 
shall consider onHs ® € N+1 are 

1 N 

H = H s ® I + I ® H R + — J2{Vi® a° + V* ® al) . (1) 

V i=l 

In [AJ1], one can find a discussion about this time renormalization and its inter- 
pretation in terms of weak coupling limit for repeated quantum interactions. 

III. The continuous limit setup 

We present here all the elements of the continuous limit result: the structure 
of the corresponding Fock space, the quantum noises, the approximation of the 
Fock space by the photon chain and [APlj's main theorem. 

III.l The continuous tensor product structure 

First, as a guide to intuition, let us make more explicit the structure of the 
photon chain. We let T9> denote the tensor product ®jn*€ N+1 with respect to 
the stabilizing sequence e . This simply means that an ortho normal basis of 71> 
is given by the family 

{e a ; a G Vjn^n} 

where 

- the set Vjn,n is the set of finite subsets 

{(ni,zi),...,(n fe ,z fc )} 
of M* x {!,..., A^} such that the n^'s are mutually different; 
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- e a denotes the vector 

Q ® . . . (g) Q <S> e it <g) ft <8> • • • <S> ft <8> e i2 <g) . . . 

where appears in ni-th copy of H, where ei 2 appears in ri2-th copy of 7i... Here 
ft plays the same role as eo in the toy model. 

This is for a vector basis on 2H>. From the point of view of operators, we 
denote by a l Jk) the natural ampliation of the operator a* to T$> which acts on the 
copy number k as a* and the identity elsewhere. That is, in terms of the basis e CT , 

a l j{k)e a = l(k,i)ea e(cr\(fc,i))u(fc,j) 
if neither i nor j is zero, and 

a l (k)e a = l(k,i)ecr e a\(k,i), 
a°j(k)e a = l(fc,o)6<T e CT u(fc,7)' 

Oo(^)e<r = l(fc,0)e«r e <r, 
where (/c, 0) £ a actually means "for any i in {1, ... , iV}, (/c, z) ^ cr". 

We now describe the structure of the continuous version of the chain of pho- 
tons. The structure we are going to present here is rather original and not much 
expanded in the literature. It is very different from the usual presentation of quan- 
tum stochastic calculus ([H-P]), but it actually constitutes a very natural language 
for our purpose: approximation of the atom field by atom chains. This approach 
is taken from [Atl]. We first start with a heuristic discussion. 

By a continuous version of the atom chain T$ we mean a Hilbert space with 
a structure which makes it the space 

$ = (g)(I7 Ar + 1 . 

R+ 

We have to give a meaning to the above notation. This could be achieved by in- 
voquing the framework of continous tensor products of Hilbert spaces (see [Gui]), 
but we prefer to give a self-contained presentation which fits better with our ap- 
proximation procedure. 

Let us make out an idea of what it should look like by mimicking, in a con- 
tinuous time version, what we have described in T$. 

The countable orthonormal basis e ff ,a G Vjn*,n is replaced by a continuous 
orthonormal basis dxo, & £ 'Pr+,Ni where Vm+ jy is the set of finite subsets of 
1R + x {1, . . . , A^}. With the same idea as for 71>, this means that each copy of 
€ N+1 is equipped with an orthonormal basis {O, dx\ , • • • , dx^} (where t is the 
parameter attached to the copy we are looking at). 

Recall the representation of an element / of T$: 

/= E /(^) e - 
ii/ii 2 = E i^)i 2 ' 
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it is replaced by an integral version of it in 

/ = / f(<r)d X „ 

J ^JR+,N 

\\f\\ 2 =f Iff da. 

This last integral needs to be explained: the measure da is a "Lebesgue measure" 
on Vpt+ n, as will be explained later. 

From now on, the notation V will denote, depending on the context, the set 
Vjn*,n or V R +^ N . 

A good basis of operators acting on $ can be obtained by mimicking the 
operators a l Ak) of We have here a set of infinitesimal operators da l j(t), i,j G 
{0, 1, . . . , A}, acting on the "t-th" copy of € N+1 by: 

dao(t) dxa = dxa dt %<£a 

da i (t) dxa = dx<ru{(t,i)} %<£o 

da l Q {t) dXa = dXa\{(t,i)} dt l( t ,i)€a 

da l j(t) dx a = dx( a \{(t,i)})u{(t,j)} l(t,i)e<r 
for all i,j G {1, . . . , A}. We shall now describe a rigourous setup for the above 
heuristic discussion. 

We recall the structure of the bosonic Fock space $ and its basic structure 
(cf [Atl] for more details and [At2] for a complete study of the theory and its 
connections with classical stochastic processes). 

Let $ = r s (L 2 (iR + , C N )) be the symmetric (or bosonic) Fock space over the 
space L 2 (M + , € N ). We shall give here a very efficient presentation of that space, 
the so-called Guichardet interpretation of the Fock space. 

Let V (= V R + N ) be the set of finite subsets {(si, ii), ■ ■ ■ , (s n , i n )} of IR + x 
{1, . . . , N} such that the Sj are two by two different. Then V = UkVk where Vk 
is the subset of V made of /c-elements subsets of M + x {1, . . . , N}. By ordering 
the iR + -part of the elements of a G Vk, the set Vk can be identified with the 
increasing simplex E fc = {0 < t x < ■ ■ ■ < t k } x {1, . . . , N} of R k x {1, . . . , N}. 
Thus Vk inherits a measured space structure from the Lebesgue measure on M k 
times the counting measure on {1, ... , N}. This also gives a measure structure on 
V if we specify that on Vo = {0} we put the measure 5$ . Elements of V are often 
denoted by a, the measure on V is denoted by da. The a-field obtained this way 
on V is denoted by T . 

We identify any element a EV with a family {a^ i G {1, . . . , A}} of (two by 
two disjoint) subsets of M + where 

ai = {s G M + ; (s,i) G a}. 

The Fock space $ is the space L 2 {V, da). An element / of $ is thus a 
measurable function / : V — > € such that 

ll/H 2 = / \f(a)\ 2 da <oc. 
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Finally, we put O to be the vacuum vector of that is, 0(a) = 5q(<j). 

One can define, in the same way, V[ a ,b] an d &[ a ,b] by replacing IR + with 
[a, b] C IR + . There is a natural isomorphism between <&[o,t] <8> $[t,+oo[ an d $ given 
by h <g) # i— > / where /(a) = /i(afl [0,i])g(ffn [i, +oo[). This is, with our notations, 
the usual exponential property of Fock spaces. Note that in the sequel we identify 
®[a,b] with a subspace of $, the subspace 

{/ G /(cr) = unless cr C [a, 6]}. 

We now define a particular family of curves in $, which is going to be of great 
importance here. Define x£ G $ by 

1 otherwise. 
Then notice that for all t G M + we have that xl belongs to $[o,t]- We actually 
have much more than that: 

Xl ~ Xl e $[ a ,t] for all s < t. 
This last property can be checked immediately from the definitions, and it is 
going to be of great importance in our construction. Also notice that xl an d X 3 S 
are orthogonal elements of $ as soon as i ^ j. One can show that, apart from 
trivialities, the curves (xl)t>o are the only ones to share these properties. 

These properties allow to define the so-called Ito integral on $. Indeed, let 
9 = {(9t)t>0i i G {1, . . . , N}} be families of elements of $ indexed by both M + 
and {1, . . . , A}, such that 

i) 1 \— > ||<7||| is measurable, for all i, 

ii) g\ G$[o,t] for all £, 

hi) Eii jo°° WdtW 2 dt < oo, 
then one says that g is Ito integrable and we define its Ito integral 

N 



±l poo 

E/ ^ 



to be the limit in $ of 



N 



EE^(4 +1 -4) ( 2 ) 

i=i jeiv 

where 5 = {tj, jElN} is a partition of M which is understood to be refining 
and to have its diameter tending to 0, and (g*)i is an Ito integrable family in $, 
such that for each i, t i— > g\ is a step process, and which converges to (<7*)j in 
L 2 (iR+ x V). 

Note that by assumption we always have that gl. belongs to $[o,t 3 ] an d Xt j+1 ~ 
X\ j belongs to ®[tj,t j+1 ]i hence the tensor product symbol in (2). 
Also note that, as an example, one can take 

9i= E r J_ r f*' +1 p *i9id8i [titti+l[ (t) 
t . es tj+i ~ tj Jt > 
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if t G [tj,tj + \], where P t denotes the orthogonal projection onto $[o,t]- 

One then obtains the following properties ([Atl], Proposition 1.4), where Vcr 
means max{s G M + ; (s, k) £ a for some k} and where cr— denotes the set a\(Va, z) 
if (Vcr, z) G a. 

Theorem 1. The Ito integral 1(g) = £V J °° g\ dx\, of an Ito integrable family 
9 = {9 x -)f=i> is the element of $ given by 

if a = 

(cr— ) otherwise. 

It satisfies the Ito isometry formula; 

(3) 



. iy /»oo 2 /»oo 

i^)ii 2 = E/ ^ = E/ M\\ dt 



In particular, consider a family / = (fi)^L 1 which belongs to L 2 (Vi) = 
L 2 (M + x {1, . . . , N}), then the family (fi(t)Q), t G iR + , z = 1, . . . , N, is clearly 
Ito integrable. Computing its Ito integral we find that 

N 



ly poo 

'(/) = £/ fimdxt 
i=i J ° 



is the element of the first particle space of the Fock space $ associated with the 
function /, that is, 



I(f)(a) = | ^ if ^ = {(M)} 



otherwise. 
Let f£L 2 (V n ), one can easily define the iterated Ito integral on 



,,(f)= [ f(a)d X<7 



by iterating the definition of the Ito integral: 

In(f)= E //•••/ fiu...,iStlT--,tn)Sld X l t \ ...d X fc. 

ii i™e{i jv}- 70 Jo Jo 

We obtain this way an element of $ which is actually the representant of / in the 
n-particle subspace of that is 

[J n (/)](cr) = / /ii,-,in(*l» •••»*»») if O- = {(*!, il)U...U(t n ,i„)} 

1 otherwise. 
Finally, for any f £ V we put 



X 



'-p 

to denote the series of iterated Ito integrals 

^ /-oo ,t„ 



iW + E E / / •••/ /ii,..,i»(*i,---»*n)ndxJi •••dxfc- 

n=l i! i„ = l' /o Jo J ° 
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We then have the following representation ([Atl], Theorem 1.7). 

Theorem 2. [Fock space chaotic representation property] - Any element f of $ 
admits a Fock space chaotic representation 

/ = / f(a) d X a (4) 
Jv 

satisfying the isometry formula 

\\ff= [ \f(a)\ 2 da. (5) 
Jv 

This representation is unique. 

The above theorem is the exact expression of the heuristics we wanted in order 
to describe the space 

R+ 

Indeed, we have, for each t G M + , a family of elementary orthonormal elements 
{O, dx\i ■ ■ ■ ■> dxt 1 } (a basis of TL) whose (tensor) products d\ a form a continuous 
basis of $ (formula (4)) and, even more, form an orthonormal continuous basis 
(formula (5)). 

III. 2 The quantum noises 

The space $ we have constructed is the natural space for defining quantum 
noises. These quantum noises are the natural, continuous-time, extensions of the 
basis operators a* (n) we met in the atom chain 

As indicated in the heuristic discussion above, we shall deal with a family of 
infinitesimal operators daUt) on $ which act on the continuous basis d\ a in the 
same way as their discrete-time counterparts a* (n) act on the e a . The integrated 
version of the above heuristic infinitesimal formulas easily gives an exact formula 
for the action of the operators a*(t) on $: 

[*?(*)/](*) = 

s<t 

K(t)f](a)= f f(o-U(s,i)) ds, 
Jo 

[«}(*)/](*)= E /((*\(«,o)u(«,j)) 

s<t 

K(t)f](a) = tf(a) 

for i, j ^ 0. 

All these operators, except a[j(t), are unbounded, but note that a good com- 
mon domain to all of them is 

P={/G$; J \a\ \f(a)\ 2 da < oo} . 
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This family of operators is characteristic and universal in a sense which is close 
to the one of the curves x\- Indeed, one can easily check that in the decomposition 
of $ ~ &[o,s] ® ®[s,t] ® 3>[t,+oo[> the operators Oj(t) — a*(s) are all of the form 

I <g> (a)(t) - a)(s))^ la t] ®I. 

This property is fundamental for the definition of the quantum stochastic integrals 
and, in the same way as for (%?), these operator families are the only ones to share 
that property (cf [Coq]). 

This property allows to consider Riemann sums: 

£ff tfc (aj(f fc+ i)-aj(f fc )) (6) 
k 

where S = {0 = to < t\ < ■ ■ ■ < tk < ■ • •} is a partition of M + , where (H t ) t>0 is a 
family of operators on $ such that 

- each H t is an operator of the form H t <g) I in the tensor product space 
$ = 3>[o,t] <E> ^[t,+oo[ ( we sa Y that H t is a t-adapted operator and that (H t ) t>0 is 
an adapted process of operators) , 

t>0 is a ste^> process, that is, it is constant on intervals: 

k 

In particular, the operator product H tk (a* (tfc+i) ~ a j(tk)) is actually a tensor 
product of operators 

H tk ®(a){t k+1 )-a){t k )). 

Thus this product is commutative and does not impose any new domain constraint 
on the operators apart from the ones attached to the operators H t and aUtk+i) — 
a l j(tk) themselves. The resulting operator associated to the Riemann sum (6) is 
denoted by 

/■OO 

T= H s da)(s). 
Jo 

One can compute the action of T on a "good" vector / of its domain and obtain 
explicit formulas which are not worth developping here (cf [Atl] for more details). 
For general operator processes (H t ) t>0 (still adapted but not a step process any- 
more) and for a general /, these explicit formulas can be extended and they are 
kept as a definition for the domain and for the action of the operator 

/>oo 

T= / H s da){s). 
Jo 

The maximal domain and the explicit action of the above operator can be de- 
scribed but also are not worth developing here (cf [A-L]). The main point with 
these quantum stochastic integrals is that, when composed, they satisfy a Ito-type 
integration by part formula. This formula can be summarized as follows, without 
taking care at all of domain constraints. Let 

/>oo poo 

T= H s da)(s), S= K s daf(s). 
Jo Jo 
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For every t G 1R + put 

<>oo 

T t = / H s l m (s)da)(s) 
Jo 

and the same for St- We then have 

/>oo /»oo />oo 

TS = / H S S S da)(s) + / T S K S da? (a) + / tf s K s daj(s), (7) 
</o Jo 



where 



_ if (i, 0^(0,0) 

^ \0 if (z,/) = (0,0). 



The last term appearing in this Ito-type formula is often summarized by saying 
that the quantum noises satisfy the formal formula: 

da l j (s)daf(s) = 6u dcij(s). 

III. 3 Embedding and approximation by the Toy Fock space 

We now describe the way the chain and its basic operators can be realized as a 
subspace of the Fock space and a projection of the quantum noises. The subspace 
associated with the atom chain is attached to the choice of some partition of M + 
in such a way that the expected properties are satisfied: 

- the associated subspaces increase when the partition refines and they con- 
stitute an approximation of $ when the diameter of the partition goes to 0, 

- the associated basic operators are restrictions of the others when the parti- 
tion increases and they constitute an approximation of the quantum noises when 
the diameter of the partition goes to 0. 

Let S = {0 = to < ti < • • ■ < t n < ■ ■ •} be a partition of M + and 
5(S) = supj|ti + i — ti\ be the diameter of S. For S fixed, define $ n = $[t„_ 1) t„], 
nG-2V*. We clearly have that $ is naturally isomorphic to the countable tensor 
product ®neiN*^n (which is again understood to be defined with respect to the 
stabilizing sequence (0) n6 ^v). 

For all neIN\ define for i,j e { 1 , . . . , N} 

e,(n) = X ); J**"- 1 E d> n , 

i / \ _ a o(^n) ~ ah(t n -l) p 
y^n tn— 1 

a){n) = Pi] o (a)(t n ) - a}(t„_i)) o P 1]; 

a o {n) = Pi]0 ^MzAh^-, 

yt-n tn— 1 

a oW = p o], 

where for % = 0, 1 and P^ is the orthogonal projection onto L 2 {Vi). The above 
definitions are understood to be valid on $ n only, the corresponding operator 
acting as the identity operator I on the others $ m 's. 
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For every a G V = Vjn*,n, define e a from the ej(n)'s in the same way as for 

e a = Q ® . . . <g> O ® e il (ni) ® O <g> . . . <g> O <g> e i2 (n 2 ) <g> . . . 
in ®neiN*'Hn- Define T$(<5>) to be the space of which are of the form 

/ = £ /(^K 

(note that the condition ||/|| 2 = J2 a ev \f( a ) I 2 < 00 is automatically satisfied). 
The space T$(«S) can be clearly and naturally identified to the spin chain T$>. 
The space T$(«S) is a closed subspace of <E». We denote by P5 the operator of 
orthogonal projection from $ onto T<E>(<S). 

The main point is that the above operators a*(n) act on T<&(S) in the same 
way as the the basic operators of T<E> (cf [API], Proposition 8). 

Proposition 3. - We have, for all i, j = 1, . . . , N 

j a l (n) ej(n) = Sijtl 

\4o = o 

f a l j(n) e k (n) = 5 ik ej(n) 
\ a) n = 

f a°(n) ej (n) = 
\a°(n)0 = e,(n) 

f ag(ra) e fc (n) = 

\ogn = o. 

Thus the action of the operators a* on the e^n) is exactly the same as the 
action of the corresponding operators on the spin chain of section II; the operators 
a l j(n) act on T$(iS) exactly in the same way as the corresponding operators do 
on T$>. We have completely embedded the toy Fock space structure into the Fock 
space. 

We are now going to see that the Fock space $ and its basic operators a*- (t) , 
i,j G {0,1,..., n} can be approached by the toy Fock spaces T$(5) and their 
basic operators a* (n). We are given a sequence (5 n ) n6 ^v of partitions which are 
getting finer and finer and whose diameter S(S n ) tends to when n tends to +00. 
Let T$(n) = T$(S n ) and P n = P Sn , for all nelN. We then have the following 
convergence result (see [API], Theorem 10), where the reader needs to recall the 
domain D introduced in section III. 2. 

Theorem 4. 

i) The orthogonal projectors P n converge strongly to the identity operator I 
on That is, any f G $ can be approached in $ by a sequence {f n ) n£]N such 
that f n G T®(n) for all n G IN. 
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ii) If S n = {0 = tft < t% < ■ • • < tl < ■ ■ ■}, then for all telR + , all ij = 
1, . . . , n the operators 

E 

k;tl<t 

E y/tk-tk-i4W, 
k-ti<t 

E ^t-n-A{k) 
k-ti<t 

and E (**-*k-iK(*) 

fc;t£<* 

converge strongly on V to a l j(t), a l Q (t), a®(t) and a[](t) respectively. 

We have fulfilled our duties: not only the space T&(S) recreates 71> and its 
basic operators as a subspace of $ and a projection of its quantum noises, but, 
when 5(S) tends to 0, this realisation constitutes an approximation of the space 
$ and of its quantum noises. 

III. 4 Quantum Langevin equations 

In this article what we call quantum Langevin equation is actually a restricted 
version of what is usually understood in the physical literature (cf [G-Z]); by this we 
mean that we study here the so-called quantum stochastic differential equations as 
defined by Hudson and Parthasarathy and heavily studied by further authors ([H- 
P], [Fag]). This type of quantum noise perturbation of the Schrodinger equation 
is exactly the type of equation which we will get as the continuous limit of our 
Hamiltonian description of repeated quantum interactions. 

Quantum stochastic differential equations are operator-valued equations on 
H s <8> $ of the form 

N 

dU t = E L)U t da)(t), 

i,3=0 

with initial condition Uq = I. The above equation has to be understood as an 
integral equation 

U t = I+ E L )U t da){t), 

the operators L l - being bounded operators on Hs alone which are ampliated to 
H s ® $. 

The main motivation and application of that kind of equation is that it gives 
an account of the interaction of the small system Hs with the bath $ in terms of 
quantum noise perturbation of a Schrodinger-like equation. Indeed, the first term 
of the equation 

dU t = L° U t dt+... 
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describes the induced dynamics on the small system, all the other terms are quan- 
tum noises terms. One of the main application of these equations is that they give 
explicit constructions of unitary dilations of semigroups of completely positive 
maps of B(Hs) (see [H-P] and also section VII of this article). 

Let us here only recall one of the main existence, uniqueness and boundedness 
theorem connected to quantum Langevin equations. The literature is huge about 
those equations; we refer to [Par] for the result we mention here. In the following, 
by coherent vectors we mean elements of the space S generated by the u <g> e(f), 
with ueH s , f e L 2 (M + ; C n ) and 

Hf)](<T)= n /*(«)' 

(s,i)ecr 

the usual coherent vectors of the Fock space 

Theorem 5.- If all the operators L l - are bounded onHs then the quantum stochas- 
tic differential equation 

U t = I+J2 / L)Usda)(s) 

admits a unique solution defined on the space of coherent vectors. 

The solution (U t ) t>0 is made of unitary operators if and only if there exist 
on Hs, a self-adjoint operator H, operators Li, i = 1,...,N and operators S*, 
i,j = 1,...,N such that the matrix (5j)i,j=i,...,jv is unitary and the coefficients 
L l j are of the form 

1 N 

L° = -(iH+-J2L* k L k ) 
k=i 




N 

L l = -Y, L l S " 

k=i 

L l j = Sj — SijI. 

m 

III. 5 Convergence theorems 

We are finally able to state the main result of [API] which shows the conver- 
gence of repeated interactions models to quantum stochastic differential equations. 

Let r be a parameter in iR + , which is thought of as representing a small time 
interval. Let U(r) be a unitary operator on Hs ®€ N+1 , with coefficients Uj(t) as 
a matrix of operators on Hs (this operator has to be though of as corresponding 
to the unitary operator U of section II). Let Vfe(r) be the associated repeated 
interaction operator: 

V k+1 {r) = U k+1 {r)V k {r) 
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with the same notation as in section II. In the following we will drop dependency 
in t and write simply U, or V/-. Besides, we denote 

1 

for all i,j in {0, ... , N}. That is, for i, j > 1 

Note that from now on we take the embeding of Tl> in $ for granted and we 
consider, without mentionning it, all the repeated quantum interactions to happen 
in H>(t), the subspace of $ associated to the partition S = {U = it; i G M}. The 
main result of [API] (Theorem 13 in this reference) is the following. 

Theorem 6.- Assume that there exist bounded operators Lp i,j e {0, . . . , n} on 
TCs such that 

UUt)-5 11 I 
lim 1— = L) 

for all i, j = 0, . . . , n. Then, for almost all t the operators V[ t / T ] converge strongly, 
when t — > 0, to V t , the unitary solution of the quantum stochastic differential 
equation 

n 

dV t = L)V t da){t) 

with initial condition Vq = I. 



IV. The G.N.S. representation of the heat bath 

In order to apply Theorem 6 to our repeated interaction model, we need the 
state of the photon to be a vector state (i.e. a pure state) instead of a density 
matrix. This is easily performed by considering the so-called G.N.S. representation 
(or cyclic representation) of the photon system. 

This representation can be described in the following way. Consider the space 
Ti = £(<D N+1 ) of endomorphisms of <U N+1 . Consider a given density matrix pp 
on Ti, which is supposed to be in diagonal form pp = diag(/3o, . . . , Pn), where all 
the Pi are strictly positive. The space H = C(€ N+1 ) is made into a Hilbert space 
when equipped with the scalar product : 

(A,B) = tr(ppA*B), 

for all A,BeH. The associated norm on H is denoted by ||-||. This Hilbert 
space is (iV + l) 2 -dimensional and we shall describe one of its orthonormal basis as 
follows. We denote by X® the identity endomorphism. Then, for i = 1, . . . , AT, we 
put X\ to be the diagonal matrices with diagonal coefficients {X], . . . , Xf} such 
that 

(Xi,Xj) = 8 ij 
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for all i,j = 0, 1, . . . N. Such a family clearly exists for its diagonal elements are 
obtained by extending the vector (1, . . . , 1) e € N+1 into an orthonormal basis of 
<D N+1 equiped with the scalar product 



N 



^ ^ Pi ^iUi 



i=0 

For i ^ j e {0, . . . , iV} we put to be the element of H given by 

It is then a straightforward computation to check that {Xj-,i,j = 0, . . . N} forms 
an orthonormal basis of H. 

We now have the usual G.N.S. representation n of C(€ N+1 ) into C{H) given 

by 

n(A)B = AB, 

for all A G /^(C^" 1 " 1 ), B EH. In the framework of that representation, note that 
Xq is then the vector state on H which represents the state pp on € N+1 : indeed, 
for all A e C(€ n+1 ) we have 

(X°,7r(A)X°)=tr( P(3 A). 

That is, in the orthonormal basis we have choosen, Xq is the only important 
vector that could not be choosen to be different. The rest of the choice for our 
orthonormal basis is just convenient for the computations, but the final result does 
not depend on it. 

Now, any operator K on Hs <8> € N+1 is transformed by 7r into an operator on 
Us ® H. That is, n(K) is a (N + l) 2 x (N + l) 2 -matrix with coefficients K$ in 
C(Hs)- These coefficients are given by 

<i = tr w (^(Xf)*Mj) 

where tr-^(H) denotes the partial trace of H along TC, that is, this is the operator 
on Hs given by the sum of the diagonal coefficients of if as a £(7ig)-valued 
(N + l) 2 x (N + l) 2 -matrix. 

When taking the continuous limit on the space Hs <E> ®in*H, we end up into 
the space 

H s ®(g)H, 

R+ 

that is, 

The associated quantum noises, following the same basis, are thus denoted by 
da% l (t),i,j,k,l = 0,...,N. 
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V. The limit quantum Langevin equation 



We are now in conditions to apply Theorem 6. We consider the repeated 
interaction model described in section II, with its associated operators H, £4, 14. 
Taking the G.N.S. representation of all that we end up in the space Hs <E> ®jn*H, 
which we embbed inside a continuous tensor product Hs <S> ®wt+H, as explained 
in section III. 3. The main result of this article is then the following. 



Theorem 7.- In the continuous limit r — > 0, the repeated interaction dynamics 
Vji/ T ] converges strongly on Hs <8> 3>, for all t, to the (unitary) solution of the 
quantum Langevin equation 



dU t = 



N N 



i=0 



i=l 



Utdt 



N 

~ i E [v 7 ^^ U t da^° (t) + JJ V* U t dal%t) 



i=i 



WPi V* U t da^t) + VPo Vi U t da^it) 



0,0/ 



(8) 



Proof 



First, we represent the operators H and U on Hs <E> € N+1 as (N + 1) x 
(N + l)-matrices with coefficients in £(Hs), following the orthonormal basis 
{e ,ei,...,ejv}. We get 



H(r) 



( H s + <y I 



1 V 2 



\ 



7^* 
Hs + iil 






l y* 


Hs + I2I 







\ 



J—V* 

N 






H s + ^nI J 



We need now to compute the associated unitary operator U = e~ irH in the 
same framework. But as we wish to apply Theorem 6, we do not need to wonder 
about the exact expression of U, but only the expansion of its coefficients in powers 
of r up to the pertinent orders given by Theorem 6. We get that U is represented 
by the matrix 



19 



( I — irHs — it^qI -iV^V* + o(r 3 / 2 ) 
-ir£f =1 ^ + (T 2 ) 

-i\fr Vi + o(r 3 / 2 ) I — irHs — ir^I 

_i T y 1 y l * + o(r 2 ) 



-iv^V^ + o(r 3 / 2 )\ 



4^ + ^ 2 ) 



V 



-^VV + o(r 3 / 2 ) -irVV^* +o(r 2 ) ... I - irH s - ir<y N I 

-iV N V*+o(r 2 )/ 

This is for the expression of £7 as an operator on Ws® € N+1 . Now our aim 
is to compute the coefficients of the matrix of tt(U) in the G.N.S. representation, 
that (N + l) 2 x (N + l) 2 -matrix with coefficients in C{Hs). 

As we already discussed in section IV, the coefficients of tv(U) are obtained 
by computing the quantities 



U% = (X?, tt(U)X]) = tr n ( Pf3 {XffUX^ 



In order to get U ' Q we have to compute 



tr H (p p X°UX°)=tr H (p p U). 



This gives, using A = 1 



u o,o 



+ itH s - ir lx I - \tV x VZ + o(t 2 )\ + . . 



u: 



0,0 
0,0 



+ f3 N \J - irH s - ir lN I - -tV n V£ + o(t 2 ) 

N 1 N 

: I - itHs ~irJ2 Pil* J2^V*V t + mV*) + o(r 



i=i 



i=0 



We now compute the U ' J terms, for i ^ j. That is, we compute 



This trace is equal to 



Ul% = te n {ppUX)). 



(ei,ppUej) = y/~Pi(ei,Uej). 



Then two distinct cases appear. If j = we have 

(I- irHs - ir 10 I - \t Z*Li ^ + o(r 2 ) \ 
-i^TVt + oiT 3 / 2 ) 



Ue = 



\ 



iVt~V n + o(t 3 / 2 ) 
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/ 



(9) 



and thus 



C^'S = -iv^>/AVi + (T 3 / 2 ). 



But when j ^ we have 



Now, if i = we get 



( -iy/fV? +o(r 3 / 2 ) 

_ i T ^* + (r 2 ) 

7 - M s - zr 7i / - \rVjV* + o(r 2 ) 

V - \tV n V* + o(r 2 ) 

^ = -^v^^+o(r 3 / 2 ), 



if % 7^ and i ^ j we get 

A similar computation gives 

C = -V^/^ + o(r 3 / 2 ), 

and, still for i ^ j and i,j^0 

U^ = -\^tV 3 V:+o{t 2 ). 

Now, consider i ^ j and /c 7^ Z, we have 

^'j = trwCp/j (Xf)*UX*) = J2( X iPP e P> UX* e p ) 

p 

1 

= 8ik{ei,Uej). 
That is, if I = and j ^ 

^; J = ^(-z^+o(r 3 / 2 )) 

and for / = j = 

JV 



C#° = 8 ik (I - itH s - ir l0 I - I £ ^ + °(^))- 

i=i 

If I 7^ and / 7^ j with j 7^ we have 

E$ = M~tW7 + o(t 2 )), 

if/=jV0 

U l k J j = 6 ik (I - itH s - ir lx I-\t V x V{ + o(t 2 )) 
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and finally if I ^ and j = then 

^;? = ^H^ + o(r 3 / 2 )). (20) 

We now study the Uj'j terms. Recall that the X\ are diagonal matrices whose 
coefficients (A°, . . . Af), i = 0, . . . , N form an orthonormal basis of (C^ -1-1 equipped 
with the scalar product £^ fiixlyi. We get 

AT 

j = tr( P/3 l/X?) = ]T(e fc , P/3 l7X*e fc ) 

fc=0 

AT 
k=0 

N 1 

= ^\ k t {I - irH s - ir lk I - —t V k V k * + o(r 2 )) 
k=i 

1 N 

+ p X°i(I ~ irH s - ir lQ I - - £ VfV, + o(r 2 )) 

i=i 

N 1 

= Pk\ k d-iri k I - -r V k V k * + o(r 2 )) 



k=i 

N 



+ /3oA°(-zr 7 o/ ~lj2v*V i + o(r 2 )) (21) 



2 

i=i 



where we have used J2 k (3 k \ h = 0. In the same way, 

N — 1 

u l? = E«(-^ J - o T w +o(r2)) 



fc=i 

1 " 



+ /3oA°(-zr 7 o/ - - £ V?^ + (r 2 )). (22) 



Finally, using £ fc=0 /3 fc A ? *i = 

8ij, we obtain 

TV 

U;'i = J2p k ^\Hek,Ue k ) 

fc=0 

* — 1 
= - zr# s ) + J^fa ^X k (-iT lk I - —t V k V k * + o(r 2 )) 



k=i 

1 N 



+ /?oAQa°(-zt 7o - + °( r2 ))- ( 23 ) 



The last terms to be considered are those of type U k \ (and conversely U i \ ), with 



k,i 
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k 7^ / and i ^ 0. We have 

Kl = trfo (XfyUXt) = J> p , P/3 ( X iTUXle t 

p 

= J2(1 p \r(Xfe p ,Ue p ) 



(3 k X^^=(e h Ue k ) 



which gives 



or 



y/^^-iy/rVi + o^ 2 )) (24) 

V^(-^t^* + o(t 2 )) (25) 

depending on k = or not. The case of U i \ is similar and needs not be explicited 
for anyway it will not contribute to the continuous limit. 

We can now apply Theorem 6. Following the rules of Theorem 6, we need to 
check that there exists bounded operators L\ 3 X on Hs such that 

Uk!i- s (i,Mk,i)I _ ij 
s - l™ 7* ~ k > 1 ' 

T ^ u f k,l 

where £q'q = 1, = £q'q = 1/2 and the others e\\ are equal to 0. 
Equality (9) shows that 

N 1 

i=l i 



By (10) we have 
and in the same way 



J 0,i 

by (11), (14) and (13) respectively. 

7-0,0 



The other terms Uq J and U i '• appear to be of order r in (12) and (15), while 



only their \fr part contributes to the limit. As a consequence Lq J = L i 'j = 0. 

Terms of the form C/^'j (equations (16) to (20)) contribute in the limit via the 

order 1 terms in U % k \ — <%,j),(fc,j) I, that is in anycase (the / term in (17) and 
(19) does not contribute as it indeed appears only when = (k, I). 

The same holds for Uq 1 and U i \ which gives Lq q = L { \ = (equations (21) 
and (22)). 

The terms U % - ,% - contribute in the limit via the order 1 terms of U\'\ — Ssjl. 

Jij JiJ ** 

Following (23) we get a null contribution in all cases. 
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Finally, equality (24) and (25) show that the last coefficients also vanish in 
the limit. 

This exactly gives the announced quantum Langevin equation. ■ 

Remark: One can only be impressed (at least that was the case of the 
authors when performing the computations) by the kind of "mathematical miracle" 
occuring here: all the (N + l) 2 terms fit perfectly in the type of conditions of 
Theorem 6. The number of cancellation one may hope for happens exactly. 

Equation (8) takes a much more useful form if one regroups correctly the 
different terms. Indeed, put 



Po 



«0 (*) + 



A- A 



A 



a o,o(*) 



Po - Pi 



Po 0,i /,\ . / Pi 0,0/,\ 



Po - P 



(26) 
(27) 



and 



W i = -iy/p -p i V i 

then the equation (8) simply writes 

N , N 

dUt = - ii 



Ufdt 



(28) 



VI. Thermal quantum noises and their properties 

In this section we concentrate on the particular quantum noises (26) and (27) 
that appeared above. We shall show that they are natural candidates for being 
qualified as "thermal quantum noises" . We also show that the form of equation 
(25) is the generic one for unitary solutions, in the thermal case. 

The situation that has appeared in the previous section can be summarized 
as follows. 

We consider the quantities /?, 7$ and thus A as being fixed. 

First of all, there is no need to consider a Fock space over L 2 (M + ; €^ N+lS> ) 
anymore, for most of the quantum noises a^ J z (t) do not play any role in equation 
(8). More economical is to consider a double Fock space: 

$ = T S (L 2 (R+, C N )) ® T S (L 2 (1R + , € N )). 

Each of the copies of the Fock space accomodates the quantum noises 

a) (t) <g> / and I <g> a) (t) 

respectively, which we shall denote more simply by 

a)(t) and b)(t) 
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respectively. 

Form the operator processes 



* = ^* ) + fiHi* (29) 



For every / G L 2 (iR + ; C n ) with coordinates (/;) in the basis {ei, . . . , e n } put 

n » 



and 



n » 



i=l jR+ 

Proposition 8.- The operators A(f),A*(g) form a non-Fock representation of 
the CCR algebra over (L 2 (M+;€ N )). 

Proof 

The operators A(f) and A*{g) have similar properties as the usual quantum 
noises. In particular, they admit a quantum stochastic integration theory, which 
is completely identical to the usual one. This does not need to be developed here. 
We shall only prove that the quantum Ito formula (7) is now driven by the rules: 



dA (t) dA° t (t) = J 30 dt 



and 

Indeed, we have 
dA (t)dA° i (t) = 



dA°i(t) dA (t) = ^ - dt. 



Po-Pi 

Pi 
Po-Pi 



<!< d4(t) + J_ft_d&J(f)) L da?(f) + , /-A- db (t) 



Po — Pi y Po — Pi 1 ) \ \l Po — Pi % V Po ~ Pi 

P ° da° (t) = dt 



Po-Pr UW PO-P 

and 

dA° i (t)dA (t) = 



! da°(t) + J db (t) ) ( J-J°— da (t) + J -J*— db°(t) 



Po-Pi lw ]j Po-Pi UW MVA)-A \l Po-Pi 

& dbl(t) = —^—dt. 



Po ~Pi Po — Pi 
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By the quantum Ito formula we get 

[A(f),A*(g)} = jr[ W)9i(t) ( §^|i) dtl 

= (f,g)i. 

In other words, the operators A(f),A*(g) form a representation of the CCR al- 
gebra over (L 2 (M + ;€ N )). But this clearly a non-Fock one for the creation and 
annihilation operator attached to this representation do not generate the whole 
creation and annihilation operators of the underlying (double) Fock space. ■ 

Now, let us form the associated Weyl operators 



We wish to compute the statistics of W(f) in the vaccum state O. For this purpose, 
we use the following notation. If H is any operator on € N+1 , then it acts on 
L 2 (R + ;€ N ) by 

N 

[Hf](s) = X + J2^f^)- 

This has to be understood as follows: in general Ao is chosen to be equal to 0, 
thus H acts on L 2 (1R + ;€ N ) as a multiplication operator. In our case it is the 
multiplication by a constant (vector). 

Theorem 9. We have 

(O, W(f) O) = exp coth(/^)/>) 

for all f e L 2 (M + ;€ N ). 

Proof 

We have 

N ' 



+ ^1^7, /,(s) i<(s) + V ^ft m w ' ' 



If we put 



and 



N roc 

(/) = £/ 



i=i ^° 
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then the above expression shows that 

A(f) + A*(f) = a(f) + b(f) 

where 



Denote by W a and the usual Weyl operators associated to the noises a and b 
respectively. We have clearly shown that 

W(f) = W a (f)®W b (f). 

As a consequence, using usual computations on the Weyl operators 

(n,w(m} = e W (-](\\f\\ 2 + f\\ 2 ^ 

|2 



= ex p (-iE( coth (/3|)ii/,n 2 ) 

= exp coth(/^)/)V 



We recover an analogue of the usual K.M.S. state statistics for a free Boson 
gas at thermal equilibrium. Let us discuss that point more precisely. Usually, 
the Hamiltonian model for a quantum heat bath is as follows. We are given a 
function u)(s) (in Fourier representation actually, by this does not matter much 
here) and the Hamiltonian of the heat bath, on the Fock space V S (L 2 (1R + ; €)) is 
the differential second quantization operator dT(u>) associated to the multiplication 
by to. 

In our discrete model, if we take the typical interacting system TH.R to be <D 2 
but with a Hamiltonian depending on the number of the copy: 

then in the continuous limit, the corresponding Hamiltonian on the Fock space 
T S (L 2 (1R + ; €)) is indeed dT(u) (under some continuity assumption on u, cf [AP2]). 

In the case we have described here, the situation is made a little more com- 
plicated by the fact that we considered a chain of € N+1 instead of (C 2 , but a lot 
easier by taking a constant Hamiltonian Hr. The time-dependent case stays to 
be explored, no doubt it will give rise to the usual free Bose gaz statistics. 
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Remark : The parameter (3 used here is supposed to be the inverse of the 
temperature of the heat bath (more exactly 1/kT). If we make the temperature 
go to 0, that is, (3 goes to +oo, then 

A) i 

Po-Pi 1 - e-/3(7i-7o) 

converges to 1, for 7^ — 70 > by hypothesis, and 

Pi e-fr* 

Po-Pi 1 - e-/3(7i-7o) 

converges to 0. This makes all the noises &*■ being useless and A(g) = a{g), for all 
g. We recover the usual quantum noises, the usual Weyl operators. This means 
that the usual quantum noises are the temperature ones. 

We now describe which kind of quantum Langevin equation, driven by those 
thermal quantum noises gives rise to a unitary evolution. 

Theorem 10. ^4 Langevin equation of the form 

n 

dU t = K°U t dt + J2 {K°U t dA$(t) + K l U t dA (t)) , U = I, 

i=0 

where the coefficients K l - are all bounded operators on Hs, always admit a unique 
solution on the set of coherent vectors. The solution is unitary if and only if it is 
of the form 

f< - (-*» - \ t {A w ' w ' + A w ' w ' )) u > dt+ 

n 
i=0 

for some bounded operators Wi, i = l,...,n, on Hs and a self-adjoint bounded 
operator H on Hs- 

Proof 

The existence and uniqueness result is a simple consequence of the one quoted 
in Theorem 5. 

For characterizing the unitarity we use algebraical (formal) computations, 
which are the same as for the proof of Theorem 5. The analytical part of the proof 
is totaly identical to the one of Theorem 5. There is no need to develop it here. 

Our equation is of the general form 

n 

dU t = K° Q U t dt + J2 i K ?U t dA^it) + Kp t dA (t)) 

i=0 

We thus also have 

n 

du; = u t *(K° y dt + J2 {u t *(K?y dA (t) + u;(K* y dA«{t)) . 

i=0 
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By the quantum Ito formula we get 

d(U*U t ) = {dU* t )U t + U* dU t + dU t * dU t 

n 

= u t *(K° yu t dt + {u t *(K?yu t dAi(t) + u t *(K* yu t dA°(t)) + 

i=0 
n 

+ U;K° Q U t dt + J2 {u;K»U t <L4?(f) + U t *K* U t dA (t)) + 



i=0 



+ TfT E u;(K°y K °u t dt + u:(K^yKp t dt 

Po - Pi ~i Po - Pi ~[ 

= u; (^K° y + k° + J2 {j^^^ + ~^zj i (^o)^o)) u t dt 

n n 

+ E ( u t ((*?)* + K o) u t ) dAi( t ) + Y {u; my + k°) u t ) dA*®. 

i=0 i=0 

By a similar computation we obtain 

d(u t un = 

= ( u t u;(K° Q y + K° u t u; + J2 (^^-^o^^(^o)*+ 

+ A K ' U ' U ' {K ' ) ")) dt 

n n 

+ Y{u t uz(K°y + Ki ) u t u;) dA (t) + + J2{u t u;(K^y + K?u t u;) dA°(t). 

i=0 i=0 

Asking both to be equal to for every t is equivalent to the following conditions : 

K = -(*?)* 

(*8r + *8 + E {irhrWrxt + jAw^^r) = °- 

~[ \Po - Pi Po — Pi J 

Put 



the last condition above exactly says 

K = -K*. 

We thus obtain the announced characterization. ■ 
The attentive reader has noticed that this is exactly the form of equation (28) ! 

VII. The Lindblad generator 

Going back to the usual Langevin equations of Theorem 5, let us recall a very 
important theorem, which is the main point in using quantum Langevin equations 
in order to dilate quantum dynamical semigroups. 
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Theorem 11. - Consider the unitary solution (Ut) t>0 of the quantum Langevin 
equation 

N N n 

dU t = -(iH+-J2 L* k L k )U t dt + J2 W t U t da^t) - ^ W£U t da^t). 



2 

k=l i=0 i=0 



For any bounded operator X on Tis, the application 

t i-> P t (X) = (O, U*(X <g> I)U t Q) 
is a semigroup of completely positive maps whose Lindblad generator is 



1 N 

C(X) = i[H, X}-- J2(W*WiX + XW*Wi - 2W*XW l ). 



i=l 



In our thermal case we have the following the form for the Lindblad generator. 

Theorem 12. - Consider the unitary solution (U t ) t>0 of the thermal quantum 
Langevin equation 

dUt = (- iH - 1 1 (^w; Wi + j^jmw;) ) v, *+ 

N 

For any bounded operator X on Hs, the application 

t h-> P t (X) = (O, Ut(X <g> I)U t Q) 
is a semigroup of completely positive maps whose Lindblad generator is 

1 N 3 

C(X) =i[H,X]--J2 tt^V ( W * W i X + XW*Wi - 2W*XWi) 

* ~[ Po — Pi 

1 N 6 

- 9 E t^V ow* + xw * w : - Wixwt) . 

1 ~[ Po - Pi 



Proof 

The basic computation is the same as for Theorem 11: 

- By the "thermal quantum Ito formula" (see proof of Proposition 8) one 
computes 

d(U*(X <g> I)U t ) = (dU*)(X <g> I)U t + U*(X <g> I){dU t ) + (dU*)(X <g> I)(dU t ); 

- The only contributing term when averaging over the vacuum state is the 
coefficient of dt, that is, we get the equation 

d(tt, U?(X <g> I)U t Q) = (O, Ut(C(X) <g> I)U t tt) dt. 

The solution is clearly a semigroup with generator C. ■ 
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Let us write down in a corollary, the Lindblad generator in perspective with 
the initial Hamiltonian. 



Corollary IS- If the repeated interaction model is having the following total 
Hamiltonian: 

1 N 

H = H s ® I + I ® H R + — {Vi ® a°i + V* <g> a l ) 

i—l 

then the associated Lindblad generator in the continuous limit is 

1 N 

C(X) = i[H s , X]--J2Po (V*V % X + XV*V % - 2V*XV i ) 

i=l 
1 N 

- g ^ x + XViV i - 2V i XV *) ■ 

i=l 



VIII. Thermalization 

In this section we answer a very natural question in this context. Consider 
a given quantum system TCs with a given Hamiltonian Hs- Is there a natural 
Lindblad generator C (in the Schrodinger picture) on TCs which admits as a unique 
invariant state, the state 

and which possesses the property of return to equilibrium for this state? By "return 
to equilibrium" we mean the following: for every inital state po? the evolution 
e tc (po) converges to the state pp in the *-weak sense, that is, 

lim tr(e tc (p )X) =tv(ppX) 

for all observable X. 

We shall prove in this section that the answer to the above question is positive, 
at least if Hs is finite-dimensional. For this purpose we recall a famous result by 
Frigerio and Veri [F- V] , in a slightly extended form due to Fagnola and Rebolledo 
[F-R]. 

Theorem 14- Let 

1 n 

C(p) = -i[H, p]~-J2 ( L *i L iP + P L * L * ~ 2L *P L *i) 
i=i 

be a Lindblad generator (in Schrodinger picture). If the commutants 

{H, Li, L*; i = 1, . . .n}' and {Li, L*; i = 1, . . . n}' 

coincide then the associated dynamics possesses the property of return to equilib- 
rium. 
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We consider Hs a iV + 1-dimensional Hilbert space, with Hamiltonian (in 
diagonal form) 

/A ... \ 
Ai ... 
H s = . . . 

\0 X N J 

Consider the Gibbs state pp = (1/Zp) e~^ Hs , it is also of diagonal form with 
diagonal elements denoted by /?o , /?i , . . . , (3n- 

We put the system Hs in repeated quantum interaction with a chain of copies 
of € N+1 with the total Hamiltonian 

N 



H = H s I + I ® H R + -L ^ (v; <g> a° + V* 

v r i= i 



4), 



where K- is the matrix 



/0 ... 1 ... 0\ 
... 



Vo o ... o/ 

with the 1 being at the i-th row and where Hr = Hs- 

Note that this means that we put the system Hr in repeated quantum inter- 
action with a chain of copies of ... itself but in the desired state. This is somehow 
very natural! 



Theorem 15 -In the continuous limit the above repeated interaction model admits 
the following Lindblab generator in the Schrddinger picture: 

1 N 

C(p) = -i[H s , p]--J2(3o (V*V lP + pV*V t - 2V iP V*) 

i=i 

1 N 



This Lindblad generator admits 



i=l 



Zp 



-pH s 



as a unique invariant state and it converges to equilibrium. 
Proof 

The announced Lindblad generator is just the dual of the Lindblad generator 
described in Corollary 13. 

Let us compute C(pp). We have 

[H s ,pp}=0. 
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On the other hand, the Vj's have been chosen so that 

Pi 
Po 

This gives 

V'* Pi rr* 

i Pa = irPbV, 
Po 

These two relation give 

V?V iP f> + ppVfVi - 2V iP V* = 2V*V iPp - 2^V l V*pp 

Pi 

and 

VWPfi + PpViV* - 2V*pV i = 2V i V*p f3 - 2§-V*V lPp . 

Po 

Hence we get the result: C(pp) = 0. 

Let us consider the von Neumann algebra generated by the operators Vi , V* , 
% = 1 ... N. It is easy to see that this is the whole B(Tls). Hence the commutant 
of this von Neumann algebra is trivial. As we always have the obvious inclusion 

{H R , V„ V*; i = 1, . . .n}' C V*;i= 1,. . .n}' 

we have equality of the two commutants and Theorem 14 applies. This gives the 
return to equilibrium property and hence the uniqueness of the invariant state. ■ 

Note the important following fact: we never used the fact that pp is a Gibbs 
state, we only used the fact that it is a function of H R . Hence the above result is 
valid for any state p which is a function of H R . 
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